Effective capacitance of a single-electron transistor 
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Starting from the Kubo formula for conductance, we calculate the frequency-dependent response 
of a single-electron transistor (SET) driven by an ac signal. Treating tunneling processes within the 
lowest order approximation, valid for a wide range of parameters, we discover a finite reactive part 
even under Coulomb blockade due to virtual processes. At low frequencies this can be described 
by an effective capacitance. This effect can be probed with microwave reflection measurements in 
radio-frequency (rf) SET provided that the capacitance of the surroundings does not completely 
mask that of the SET. 
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A single-electron transistor (SET), shown schemati- 
cally in Fig. [TJ is one of the most widely studied compo- 
nents of nanoelectronics today. Numerous applications 
include charge detection in mesoscopic structures [l|, 
thermometry and single-electron pumping Q. Due 
to its high sensitivity to charging effects, a SET is an 
ideal structure to study and characterize single-electron 
and quantum effects such as Coulomb blockade and tun- 
neling. 

The orthodox theory developed by Averin and 
Likharev 0] has been extremely successful in describing 
the dc properties of a SET. More recently, however, SET 
has also been used at finite frequencies as an accurate 
electrometer in the form of rf-SET []]]. In this applica- 
tion, the impedance of the SET has been assumed to be 
completely resistive. Previously the response of a SET 
to an ac field has been studied in the framework of scat- 
tering matrix formalism 0, Q and with the Tien-Gordon 
approach @, B 0|- In this Letter we utilize the Kubo 
linear response formula to calculate the finite-frequency 
admittance of a SET in the presence of Coulomb inter- 
action. In addition to the familiar sequential tunneling 
effects, at finite frequencies the lowest order approxima- 
tion describes some virtual tunneling processes as well. 
These give rise to a nonzero reactive response even under 
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FIG. 1: (color online) Schematic drawing of a single-electron 
transistor considered in this Letter. Left electrode is con- 
nected to an ac voltage source of angular frequency u. Cen- 
tral island is coupled capacitively to a gate electrode, and the 
gate charge can be adjusted with the gate voltage Va- 



Coulomb blockade. At low frequencies this contribution 
is described by an effective capacitance, which can be 
tuned with the gate voltage. We also discuss the possibil- 
ity of detecting this gate-dependent capacitance through 
the tuning of the resonant frequency in an LC circuit. 

Usually it is the familiar geometric capacitance that 
dominates the total capacitance of a given system. There 
are, however, additional contributions that have a micro- 
scopic origin. These are the density of states capacitance 
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111 . Il2j , which results from the finite amount of ki- 



netic energy that is required to introduce an additional 
electron to a conductor, and the correlation capacitance, 
which comes from the correlated motion of electrons. 

The Hamiltonian for the single-electron transistor is 
H = Hl + Hj + Hf> + H Tl where H L = J2 V ^4 c f de- 
scribes the noninteracting electrons with eigenstates v in 
the left lead and Hr is a similar expression for electrons 
in the right lead. The island is described with the Hamil- 
tonian 



+E(N), E(N) = E c N 2 -eV G N, (1) 



where N is the number of excess electrons on the island 
and Eq = e 2 /2C is the charging energy for an island of 
capacitance C — Cl+Cr+Cg- For the sequel, it is useful 
to define 5E% = E(N ± 1) - E(N) = (±2N + 1)E C T 
eVc- The tunneling Hamiltonian, describing the charge 
transfer processes between the island and the leads, is of 
the usual form 



H T —H T l + h tl + Htr + H TR , 



H Ta = (H+ a y, (2) 

up, 

where t ufl is the tunneling matrix element between two 
corresponding states and a € L, R. The Kubo formula 
for conductance reads [TJl 



(3) 



where the retarded current-current correlation function 
is the Fourier transform of 



C? I (t-t') = -i6(t-t')([I(t),I(t')]) 



(4) 
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and the second term involving positive constant Co is the 
so-called diamagnetic term, which cancels the zero fre- 
quency divergence in the imaginary part of conductance 

The particle current operator for a single-electron tran- 
sistor can be found as the time derivative of the particle 
number in, say, the left lead. The Heisenberg equation 
of motion for the number operator yields 

N L (t) = i [H, N L (t)} = i[H T , N L (t)} , (5) 

which leads to 

I(t) = N L (t) = -i (H+ L (t) - H TL (t)) . (6) 

The current-current correlator is then 

C?j(t -t') = 2 Im 6{t - < [H TL (t), H+ L (t')] ) , (7) 

where the time-evolution of the operators is determined 
by the full Hamiltonian and the expectation value should 
be calculated in the presence of tunneling. Transforming 
this expression to the interaction picture and expanding 
to the lowest non-vanishing order in t v)1 we find 

C#(* - =2 Im 9(t - E E £mW 

x([4Wc„W,4,(Ov(t')])o, (8) 

where the subscript zero means that the expectation 
value and time-evolution of the operators should be eval- 
uated with respect to the Hamiltonian Hl + Hi + Hr. 
Using standard methods of many-body theory, the ex- 
pectation value in Eq. JHJl can be decomposed to Fermi 
functions and time-dependent exponentials with the help 
of the finite-temperature Wick's Theorem. Thus, in fre- 
quency domain Eq. (jSJ takes the form 





FIG. 2: (color online) Real (a) and imaginary (b) part of the 
zero-temperature admittance as a function of the frequency of 
the driving signal and the gate voltage. If the SET is assumed 
to stay in the ground state at all times, the patterns repeat 
periodically as functions of the gate voltage. 
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uj — e u + Efj, + 5E^ uj + e v - e M — 5E^ + irj 

(9) 



where 77 is a positive infinitesimal quantity. We transform 
the sum to an integral and using l/(x + irj) — V 1/x — 
iir5(x) the imaginary part takes the form 
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and the real part, without the diverging term cancelled 
by the diamagnetic term in Eq. @, is given at T — by 
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(11) 



Here Rt — h(2ire 2 \tv f _ l \d,Ldn)~ 1 is the usual tunneling 
resistance determined by the tunneling amplitude and 
densities of states d,L,dn at the Fermi level. Numeri- 
cal integration of Eq. |9]) shows that the temperature 
dependence of the real part is exponentially weak when 
ksT <C Ec- The conductance follows from Eq. |3]). Note 
that the real part of the conductance is obtained from the 
imaginary part of the current-current correlator and vice 
versa. 

Real part of the admittance, shown in Fig. [Ha), ex- 
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hibits the familiar Coulomb blockade as can also be seen 
from Eq. (|10[) . If the frequency of the driving signal is not 
sufficiently high to provide the required charging energy, 
real part of the admittance is exponentially suppressed 
and vanishes completely at zero temperature. For higher 
frequencies the blockade is lifted, which corresponds to 
photon assisted tunneling. When eVc — (2N ± l)Ec, 
two lowest charge states are degenerate, and the charg- 
ing energy vanishes. Imaginary part of the admittance 
is shown in Fig. EJb) . It is nonzero even under Coulomb 
blockade because of the possibility of electrons to tunnel 
back and forth to a virtual state on the island. It should 
be noted that for dc-response it is necessary to take the 
second order approximation in Ht for virtual processes to 
appear [l5| . The magnitude of the imaginary part grows 
linearly at low frequencies, as can be seen by expanding 
Eq. (fTTj) near co = 0. We obtain a linear admittance 
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which implies capacitive behavior with an effective ca- 
pacitance 



\2ttJ \R t 



2E C 



\{2N + 1)E C - eV G ][(l - 2N)E C + eV G ]J ' 

(13) 
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FIG. 3: (color online) The total capacitance of the SET 
as a function of Rt/Rq- The different curves correspond 
to different values of gate voltage (from top to bottom): 
eVa/Ec = (blue), eV G /E c = 0.2 (cyan), eVa/E c = 0.4 
(green), eVa/Ec = 0.6 (magenta) and eVa/Ec = 0.8 (dark 
green). The dashed red line corresponds to the geometric ca- 
pacitance without the quantum correction. The first order 
approximation is no longer valid when Rt/Rq approaches 
unity. Inset shows the gate dependence of the total capaci- 
tance for Rt/Rq = 2. 



where Rq — h/e 2 is the quantum resistance. Capaci- 
tance can be tuned with the gate voltage and it diverges 
at the charge degeneracy point. This divergence, appear- 
ing only at zero frequency, is an artifact of the first order 
approximation and is cancelled by higher order contribu- 
tions. Similar result was found for a single tunnel iunc- 
tion under dynamical Coulomb blockade in Ref. [lq . 

The current operator that we have used in the calcula- 
tions takes only into account the tunneling of individual 
electrons. In problems involving a time-dependent driv- 
ing it is also crucial to take into account the displace- 
ment current originating from the continuous displace- 
ment of electronic charge to satisfy the current conserva- 
tion H, [l7[ . We consider the following extension to our 
model: In addition to a "particle current channel" formed 
from the tunnel junctions we have a parallel "displace- 
ment current channel" formed from the capacitances of 
the tunnel junctions. We assume a left-right symmetric 
SET geometry with Co -C Cl, Cr. Because the gate 
capacitance is negligible compared to the other two ca- 
pacitances, the current in the left lead equals the current 
in the right lead 17f ■ Thus the total admittance is given 



by a sum of the geometric and the tunneling contribution 
G = Gq — iujC '/4, where Gq is the admittance calculated 
from the Kubo formula. The quantum correction is usu- 
ally smaller than the geometric contribution, but its gate 



dependence can still be observed provided that the tun- 
neling resistance is not too large. The total capacitance 
of the SET, Ctot = C/4 + C, as a function of R t /Rq is 
shown in Fig. [3] 

The reactive impedance can be utilized as an electrom- 
eter in a following setup (shown in Fig. [2J , similar to the 
rf-SET scheme: A resonator circuit formed of an inductor 
of inductance L, the SET and a stray capacitance Cm is 
fed an rf-signal through a transmission line of impedance 
Zq. The gate charge can then be probed by measuring 
the phase of the reflected signal. The phase of the reflec- 
tion coefficient r = (Z — Zq)/(Z + Zq), where Z is the 
total impedance of the resonator circuit, as a function 
of frequency and gate voltage is shown in Fig. [51 Phase 
changes sign at the resonance frequency, which can be 
tuned by the gate voltage. Note that in typical rf-SET 
measurements so far the parameter range has been signif- 
icantly different from those where this effect is observed, 
and thus it has stayed undetected. 

Our analysis is limited by three factors. First of all, 
the exact charge state of the island is not calculated self- 
consistently, rather, it is assumed that the island is at 
all times in the ground state. This assumption should 
be valid at least when |7iu;| < 5E^,SE^, which is satis- 
fied for most practical applications. Second, non-linear 
response effects are neglected within the applied linear 
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FIG. 4: Proposed measurement setup for the measurement of 
gate charge with the reactive impedance of a SET. 



dependent SET admittance starting from the Kubo for- 
mula. We found novel features, most notably that the 
reactive part is nonzero under Coulomb blockade due to 
virtual processes and can be described by an effective 
capacitance at low frequencies. Our results can be ex- 
perimentally demonstrated by rf-reflection measurements 
which should reveal a strong gate voltage and frequency 
dependence. 
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P. J. Hakonen and E. B. Sonin for discussions. TTH is 
supported by the Academy of Finland. 
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FIG. 5: (color online) Phase of the reflection coefficient in 
a setup shown in Fig. [4] Phase changes sign at the reso- 
nance frequency, which depends on the gate voltage. The 
parameters that were used are Rt/Rq ~ 2, EcL/TiRq — 0.6, 
C||/C = 250 and Z /R Q = 0.002. For E c /k B = 1 K, these 
correspond to Z ps 50 ft, L » 160 nH and C|| » 250 fF. 
Same values were used in the measurement of Ref. [To|| . The 
resonance frequency at zero gate is around ljo ~ 6 x 10 9 s _1 . 

response theory. Third, near charge degeneracy points 
and for Rt/Rq ;$ 1 second and higher order effects be- 
come relevant. However, none of the above mentioned 
limitations prevent the possibility to experimentally ver- 
ify our analysis with currently available technology. 
In conclusion, we have calculated a frequency- 
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